ON THE UNIFORM CYCLE-ROOTED SPANNING TREE IN 1? 
ADRIEN KASSEL, RICHARD KENYON, AND WEI WU 

Abstract. We compute the asymptotics of the first and second moments of the 
fvj area of the cycle of a random cycle-rooted spanning tree (spanning unicycle) of any 

^-H sequence of graphs Qn C Z^, such that -Qn approximates a bounded domain £) C 

^^ C We show that the first and second moments grow like - log n and C- Area(D)n^, 

respectively, for an explicit constant C — C{D). 

We use these results to give a lower bound for the first and third moments of 

the length of the random loop obtained by adding an independent random edge to 

a uniform spanning tree on C/„. 
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1. Introduction 

A cycle-rooted spanning tree (CRST), also called spanning unicycle, on a graph 
^ is a connected spanning subgraph containing a unique cycle (in other words a 
connected subgraph with as many edges as vertices). More generally a cycle-rooted 
spanning forest (CRSF) is a subgraph each of whose connected components has a 
unique cycle, that is, is a cycle-rooted tree. 

On any connected finite planar graph Q we denote by Tg the random variable 
corresponding to the uniform measure on CRSTs. We let Lg and Ag be the com- 
binatorial length (number of edges) and the combinatorial area (number of faces) of 
the cycle of Tg. 
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Figure 1. A sample of a uniform CRST on a 21 x 21 square grid 
conditioned on having a longish cycle. 



We show (Theorems p^ and H| that for the n x n square grid ^ the area satisfies 
E(Ag) = ^ log(n) + o(log(n)) and E(A2) = Cn^ + o{n'^) where C = ^-^ ^ .281 and 
(5 is the constant 



/? 



E 



fc,Z=l, both odd 



m^{k^ + p 



.527, 



We extend this result (Theorem ^ to any exhaustion Qn of Z^ such that ^Qn 
approximates a bounded domain D G C We show that the expected area satisfies 
the same asymptotics, that is E(Ag^) = ^\ogn + o(logn), and that E(AgJ grows 
hke C{D)\D\n'^, where l^l is the area of D and 



depends on D as the integral of its Dirichlet Green's function g%. We call C{D) the 
mean normalized exit time of D, since it is (when multiplied by |-D|/4) the mean 
exit time for Brownian motion started from a uniform random point; it is invariant 
under similarities of D. For an n x m rectangle with aspect ratio m/n — )■ r, the 
mean normalized exit time is C = — ^^r-^, where 

7T-0 1 



/3(-) = E 



k,l=l, both odd ^ ^ 

For the disk the mean normalized exit time is C = l/vr, see below. This is the 
maximal possible value of C{D) among all Jordan domains. The constant C{D) also 
has a mechanical interpretation as C{D) = t^P{D), where P{D) is the torsional 
rigidity of a cylinder beam with cross-section Z), see Section |5j 

The uniform CRST is closely related to the uniform spanning tree, but on general 
graphs no formula is known even for the total number \g of CRSTs. On increasing 
exhaustions of the square grid {Qn)n>i the number of CRSTs was recently shown by 
Levine and Peres ^2J (based on results of Kenyon and Wilson [D]) to be ^-^Kg^{l + 
0(1)), where \Vn\ is the number of vertices of Qn and Kg^ is the number of spanning 
trees of Qn- As a simple corollary they showed that the expected length of the 
cycle in a uniform CRST on the grid is 8 + o(l). A similar argument shows that 
E(Lg^) = 8 + 0(1), when -Qn approximates any bounded domain D C C. 

The unusual feature of this random cycle, that the average length is finite while 
the average area grows, is a consequence of the "long tail" of the distribution. It is 
natural to conjecture that the probability that the area is K decays like K~^, and 
(based on the known relation between length and diameter for the branches of the 
uniform spanning tree [7]) the length L and area K are related by L = 0{K^/^). 

Using the techniques of [Ij, one can explicitly compute the probability that the 
cycle has any given shape (in particular any given length), in the limit n — )■ 00, since 
removing a fixed edge from the cycle leaves a uniform spanning tree containing the 
remaining set of edges. The probability of lengths 4, 6, 8 are respectively 
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2. Laplacian, Green's function, and transfer impedance 

Let Q he a finite graph with vertex set V. We let Kg be the number of spanning 
trees and Xg the number of CRSTs of Q. For any edge e (respectively any face /), 
we write e E Tg (respectively / G Fg) to express the event that the edge is part of 
the cycle of Fg (respectively, that / is in the interior of the cycle of Fg). 

The main tool for computing the quantities we are interested in will be the transfer 
impedance. Let us recall some basic facts about it and its link with the Green's 
function. 

We define the Laplacian A : M^ ^ M^ by Af{v) = 'Zv'^v /(^) ' /(^'). the sum 
being over nearest neighbors f ' of i;. As long as Q is connected A has kernel consisting 
of the constant functions. A is invertible on the orthocomplement Vq of the constant 
functions and we define the (Neumann) Green's function G to be this inverse. By 
abuse of notation for vertices a,b,x we denote G{a,x) — G{b,x) the function of x 
defined by A^^{6a — 6b)- 

Let A|j'j be the determinant of the submatrix of the Laplacian in which rows a, b 
and columns c, d have been deleted. 

Lemma 1 (See e.g. [T3]). 

(1) i-ir^'+'+^A'^^, = Kg (G(a, c) - Gib, c) + Gib, d) - G(a, d)) . 

Given two oriented edges ei = ab and 62 = cd, define the transfer impedance 
Tg{ei, 62) by T^ab, cd) = Gia, c) —Gib, c) — Gia,d) + Gib, d). The transfer impedance 
has a random walk interpretation and an equivalent electrical interpretation, see [T] , 
El, or 



Lemma 2 (See e.g. [13]). If we interpret the graph as a resistor network with unit 
resistances on the edges, then T{ab, cd) is the amount of current crossing edge cd 
when one unit of current enters at a and leaves at b. 

One can give an expression for T involving an eigenbasis of the Laplacian. We 
denote by (■, ■) the usual scalar product on Vq, coming from restriction of the stan- 
dard scalar product on M^. Let {fk} be an orthonormal basis of Vq consisting of 
eigenvectors for A associated to eigenvalues A^ 7^ 0. Then, for any two functions 
/,5' e Vo, we have 

/f A-i^\ _ \^ {fk,f){fk,g) 

k ^'^ 
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since this equality holds when / and g are replaced by any elements of the basis 

{fk}. 

Therefore we may write the transfer impedance as 

T{ab,cd) = {6,-6d,A-\6a-6,)) 

_ sr^ {fk, 5c - Sd) {fk, Sg - Sb) 

~ r ^^ 

ifkic) - fk{d)){fk{a) - fk{b)) 



(2) = E 



k ^^ 

Recall also the Matrix- Tree theorem 

Theorem 1 ([lU])- For any row a and column h, Kg = (— 1)''^^A^. 

3. Estimates for general finite graphs 

3.1. Moments. The general strategy for computing moments of Ag is to use the 
following simple observation. 

Lemma 3. For any integer k > 1, 

E(A^)= J2 nfu--- Jk^Tg) 

/l>--- Jk 

where the sum is over all ordered k-tuples of bounded faces. 

In this section, we use the line-bundle Laplacian j8] to derive the probabilities 
I^l/i)''' , fk £ Tg) for fc = 1 and 2. Let us briefly recall the definition of this 
Laplacian and the result relating its determinant to the CRSFs of the graph. 

3.2. Line-bundle Laplacian. Let Q = {V,E) be a finite graph. A C-bundle is a 
copy Cy of C associated to each vertex v E V. The total space of the bundle is the 
direct sum W = 0^gy C^,. A connection $ on VT is the data consisting of, for each 
edge e = vv', a complex linear isomorphism ip^yi : C^, — ;■ C^i. That is, we assign 
to each oriented edge e a non-zero complex number (fe such that y^-e = fe^- ^^ 
say that ipe is the parallel transport of the connection over edge e. We say that two 
connections $, $' are gauge equivalent if there exist non-zero complex numbers ipy 
such that ipy'ipyy' = ip'^^/ipy, that is, $' is obtained from $ by changing the basis of 
the C^,. 

The monodromy of a connection around an oriented cycle Ui, • • • , Vn+i = vi is the 
complex number nr=i fv,+uv, e C*. 

5 



We let A = A$ be the associated Laplacian acting on f G W defined by, for eacfi 
vertex v 

A{f){v) = J2fiv)-^v'J{v'), 

wfiere tlie sum is over neiglibors v' of v. 
Theorem 2 (|ll|S]). We have 

det(A)= J2 n(2-^-lM, 

CRSFs cycles 

where u is the monodromy of the connection around the cycle, for either choice of 
its orientation. 

3.3. Two technical lemmas. We first start with two technical lemmas. 

Lemma 4. Let M he a square matrix and Z a block- diagonal matrix of the same 
size: 

(Bi \ 

Bo 



B, 



\ 



0/ 



V. 



consisting of blocks of the form Bi = \^ ^r\] ' "^here {pi, qi}i<i<k are 2k variables, 

and elsewhere. We denote by Xi,yi the two consecutive rows (and columns) indices 
corresponding to the i-th block. Let F := Fm{pi, ■ ■ ■ ,Pk,qi,- ■ ■ , Qk) be the polynomial 
det(M + Z) E 'C[{pi,qi}i<i<k]. We denote by H the orthogonal projection of F on 
the space of polynomials of total degree at most 2. Then 



H = Y.i'g^My:-p,M^:-p,q,M::;^:) 



4 = 1 



+ 



Z^ 



{p^PjM^::^; + q^q,Myii:i^ + p^q,M^:^l + q^P,M^^D 



where M^ denotes the determinant of the submatrix of M where the rows indexed by 
A and the columns indexed by B are removed. 

Proof. We compute explicitly the coefficient of each of the monomials that can appear 
in H: pi, q^, piPj, qiqj, p^qi, p^q^, and qiPj, for i, j = 1, ■ ■ ■ ,k and i < j. 
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When expanding the determinant along the yi-th row, the term pi appears only 
once, and in order to get the degree 1 term, we set pj, Qj to be zero in the correspond- 
ing minor. This yields that the coefficient of Pi in F is -^[[pj]] = — M^'. Similarly, 
expanding along the yi-th column, we obtain that -F[[q'j]] = —M^\ 

We also obtain that for all i, F[[piqi]] = — M^^*. And for i < j, we have -^[[piPj]] = 
M,7,S^ and F[[q,q,]] = Mjf;*:, and F[[p,q,]] ='<;*, and F[[q,p,]] = Mi^;,f . The 
result follows by summing all the terms. D 

The next lemma is a generalization of the previous one to a degenerate setting. 

Lemma 5. Let M be a square matrix and {pi,qi}i<i<k be 2k parameters. Let Z be 

/O p. 0~ 

a block-diagonal matrix of the same size as M with blocks of the form \ qi Pi 

\0 q^ 

for i < r, and I * J for i > r + 1, and elsewhere. The corresponding rows 

and columns are labelled by Xi,yi, Zi for 1 < i < r and Xi, yi for r -\- 1 < i < k. Let 
F := det(M + Z) and H be its part of total degree at most 2. Then 

k 
i=l,{ai,bi) 

+ E (p'P^KX + ^^^X:;^ + p^<lJK::a: + mM^Si) 

^^'i-<j^k,{ai,bi),{aj ,bj) 

r 

i=l 

where there is an extra summation over {ai, hi) for each i < r for which {ai, hi) = 
{xi,yi) or {yuZi); if i > r, {ai,bi) just means {xi,yi). 

Proof. The proof follows by expanding along rows and columns for each rows/columns 
indexed hj Xi,yi,Zi. D 

3.4. Probability estimates. We won't need the following lemma to compute mo- 
ments, but state it for its relevance to the quantities computed and its own interest. 

Lemma 6. Let Q be any connected finite graph. The probability that an edge e is in 
the cycle of Vg is 

P(eGr,) = ^(l-T,(e,e)), 
Xg 

where Tg is the transfer impedance. 
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Proof. Consider the bundle Laplacian with parallel transport q ^ 1 on edge e (with 
either orientation), and 1 on all other edges. We have det Ae(g) = (2 — g — l/q)Ne, 
where iVg is the number of CRSTs containing e in its cycle. Hence 

WeGT ) - ^^- ^ detAe(g) 
Ag Xg{2-q- l/q) 

= y^(l-T(e,e)), 
Ag 

using Lemma 111 (note that Ap differs from the standard Laplacian by a single 2x2 
«oc. (, _V ^ > --- 0^ -^ "- --'-- '— 

The next two lemmas give the probabilities P(/ G Tg) and P(/i, /2 G Fg) in terms 
of the potential theory properties of the graph. 

Lemma 7. Let Q he any connected planar finite graph. Let f he a hounded face. Let 
7/ he a simple dual path from f to the unhounded face. Denote hy E{'yf) the set of 
edges the path 7j crosses, oriented counter-clockwise around f . Then 

(3) P(/Gr,) = ^ E (l- E Tgia,b)) = 'j^ Yl Tg{a,b) , 

where 7/ is a dual path from f to the boundary such that for each edge e = xy in 
E{'~ff), the path 7^ U 7/ disconnects x from y in Q. 

Proof. For a face f oi Q, and a real or complex parameter g 7^ 0, define a connection 
on a discrete line bundle over Q by taking the parallel transport equal to the identity 
on all edges except those crossed by 7^; on these, take parallel transport q (oriented 
consistently in such a way that the monodromy on any counter-clockwise loop around 
/ is g). See Figure [2| Let Aj(g) be the corresponding bundle Laplacian. 

We have 

00 

detAf{q) = Y,Nk{2-q-l/q)', 
fc=i 
where Nk is the number of CRSFs with k cycles each of which winds around /. Thus, 
the probability of the face / being inside the cycle of a uniform CRST is 

... . ^ N ^1 1 , det A f(q) 

(4) p(/erg) = -^ = — hm ^^^^ 



Xg Xg q^l2- q- l/q 

Let Cj = Xiyi be the elements of the set E{jf) and k = |£'(7j)| its cardinality. Let 
A be the matrix of the combinatorial Laplacian where we have chosen to order the 
vertices of Q starting with xi,yi,X2,y2,- ' ' i^kiVk- Let pj, qi be 2k variables and let Z 




Figure 2. Zipper connection with parameter q for a fixed face /. 



be the block-diaeronal matrix, whose first k blocks are I J: | and zero elsewhere. 

Now set F to be the determinant of A + Z. This is a polynomial in the variables pj, 
qi. Let H be its degree-at-most-two part as in Lemma |4} 
By Lemma |4| 



if = 5^(-g.A^:-p.A--p,g,A-*) 



1=1 



E (m A^::^ + q^qj^^^ + P.giA^;;^^ + g.p^A^-^^) 



l<i<j<fc 



Using the symmetry of the matrix A, this yields 



^ = E(-(p^ + ^^)^^:-p^^^^"*) 



4 = 1 



l<i<j<fc 

Let e > be small, and set pj = 1 — g = — e, q^ = I — l/q = e — e'^ + C){e^). We 
now view F and iif as functions of e. 

Recalling Theorem [T| F{e) = H{e) + 0{e^), where 

k 

l<i<j<k 
9 



i=l 



Using Lemma [T] and the definition of T, 



and 



K-g{Gix^,yj) + G{xj,y^) -G{xi,Xj) -G{y^,yj)) = -KgT{ei,ej) 



J J ' 



^Zfi = ^G {Gix^, Xi) + G{yi, yi) - G{xi, yt) - G{xi, yi)) = ngT^a, d) . 
We obtain tliat 

k k 

det A^(l - s)/{-e^) = /^g E (l - E ^(^- ^^•)) + ^(^) ' 

Tlie first equality in ^ now follows from Q. The second equahty follows from 
Lemma [2l D 

Lemma 8. Let Q he any connected planar finite graph. Consider two bounded faces 
f and f and two non-intersecting simple dual paths 7/ (resp. 7/'j from f to the 
unbounded face (resp. from f to the unbounded face). Let E{'~ff),E{'~f'p) be the sets 
of edges crossed by 7^ and 7/', oriented counter-clockwise around f,f' respectively. 
Then 



/'ere 



Kg 

a7 



Yl Tg{a,b). 



aeE('Yf),beE{Tf,) 



Proof. Take a connection with monodromy q around / and q' around /' as before, 
supported on E{^f) and E{^fi). See Figure Isl 




Figure 3. Zipper connection with parameters q and q' for two fixed 
faces / and /', where the paths potentially have common adjacent 
vertices. 
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We have 

det Afjiq, q') = ^ iV,,,,,(2 - q - l/g)^(2 - q' - l/q')\2 - qq' - l/{qq')r , 

a,b,c 

where A^a,6,c is the number of CRSFs of Qn with a + b + c cycles, a of which contain 
/ but not /', b of which contain /' but not /, and c of which contain both / and /'. 
Therefore, P(/, /' G r„) = iVo,o,i/Ag. 
When g — !■ 1 we have 

det Afjiq, q) = (iVi,o,o + iVo,i,o)(2 - q - l/q) + iVo,o,i(2 - q^ - l/q^) + 0((g - 1)^) . 

As g — )■ 1, we get 

hm ^f^f^f'^^f = N,,, + iVo,,o + 4iVo,o, . 
q^i 2 — q — l/q 

Using the previous section, we also get that 

,. detAff/(g, 1) 

g^i 2 - q - l/q 

and 

detA;,;,(l,g) _ ^ , ^ 
g^i 2 — q — l/q 
Hence 

iVoo, = Ifhm d^^^/.H^-;) - hm det A/,Hg,l) _ 1-^ detA.Alq) 
"'"•' 2Vg^i 2-q-l/q q^i 2-q-l/q 9-^1 2-q-l/q 

Let us first deal with the case where the paths 7/ and 7// don't have any common 
adjacent vertex. In that case, we deduce from Lemma [7] that the contribution is 
— Kg/Ag ^- . T(ej, e^), where the sum is over all Cj and e'j in E{'~^f) and E{jfi), 
respectively. 

Let us now consider the case where the paths 7/ and jf are adjacent to common 
vertices, say r of them. We label these r 2-edges-paths by ordered triplets of vertices 
Xi, Hi, Zi for i < r. The rest of the edges of E{jf) are labelled Xiiji for i e / and the 
rest of the edges of -E(7/') are labelled x'j, y'j for j & I'. 

The Laplacian with parallel transport q on E{pff) and q' on E{^fi) can be written 
as A + Z where Z is the following sparse block-diagonal matrix. The non-zero blocks 

/O Pi 0\ 
of Z consist in r blocks of the form \ qi Pi\ indexed by rows/columns Xj, j/j, Zi 

VO q^ 0/ 

for i <r and A; — r blocks of the form I * 1 for r + 1 < i < A;, where we ultimately 

replace pihj 1 — q = —e and g^ by 1 — l/q = e + 0{e^). 
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Using Lemma [5} we can compute the degree at most 2 part of the differences of 
four determinants of Laplacians det(A + Z) seen as polynomials in the variables Pi,qi. 
After cancellations, the remaining terms of degree at most 2 are for pf, qf for i < r; 
PiPjililjjPilj cLnd QiPj, for i < j and appropriate pairings. 

We replace Pi = 1 — q = —e and g^ = 1 — 1/q = e + O(e^). The total contribution 
(dividing by 2 — g — 1/q = —e^) is 

where the right hand side was obtained using Lemma [T] and the definition of T. 
Dividing by 2 we get the result. D 

Remark 1. It is worth noting that the right-hand sides of Lemma\2\ o,nd Lemma\^ 
do not depend on the particular choice of paths 7/ and -yfi . 

4. Estimates for finite square regions in I? 

Let Qn be the n x n square grid graph. Its vertices are indexed by [l,n] x [l,n] 
and edges join nearest neighbors. We denote by Kn and A^ the number of spanning 
trees and cycle-rooted trees of Qn, and by r„, A^ and L„ the uniform CRST on Qn, 
its area and its length. 

4.1. Explicit expression for the transfer impedance. Let Vq denote the space 
of zero-mean real- valued functions on the vertices of Qn- 
The family of functions 

2 

fkA^, y) = -JT—-. ^i V cos(7rA;(x - l/2)/n) cos{7rl{y - 1/2) /n) 

n[l + lfc=o + l«=o) 

for (A;, /) G {0, ■ ■ ■ , n — 1}^ \ {(0, 0)} forms an orthonormal basis of Vq which con- 
sists of eigenvectors for A|yg corresponding to the eigenvalues A^ ; = 4(sin^(|^) -|- 

Hence, using ([2]), we obtain that for any two horizontal edges Cj = ((x, t/j), (x + 
l.Vi)) and e^. = ((x',^^), {x' + l,yj)), 

,^ 1^ . (^_kx^ ■ . 7rfca:\ sin^(i^) cos(^^^fa^) cos{^^^^^^) 
^(e-S)-^2 2^sin( ^ )sm( ^ ) gi^2(^^/(2n)) + sin2(7r//(2n)) 

^^) n-l 

2 v-^ ,Txkx , ,Txkx\ 



-F ^j > sin( )sin( ). 

fc=i 

In particular, if x = x', we have 



12 



(6) r(e.e-)-l + lV.in'(^) "°'(i^^-(^^'^^)-(^^'^^) 



since for any x 7^ mod n, we have X]fc=i ^^^ 






4.2. Expected area. Recall that the expected area of r„ is 

E(A„) = 5^p(/Gr„), 
/ 

where the sum is over all {n — 1)^ bounded faces of Qn- 

The probabilities that appear in the sum may be computed with the help of the 
following lemma which is a special case of Lemma [7| 

Lemma 9. Let f be a bounded face of Qn one? denote by k the y-coordinate of its 
lower vertices. Denote by Ci = aibi the k east-oriented edges below f , enumerated 
from bottom to top. Then 

k 

P(/er.) = g(A:-^T(e.,e,)), 
Equivalently, 

k 



P(/Gr„) = ^5^n, 



An . , 
1=1 



where Uf{i) is the amount of current that flows left of f when Qn is considered as a 
resistor network with all edges having resistance one, and when one unit of current 
flows into Oj and out of hi. 

Proof. We consider a vertical zipper of edges below face /, see Figure |4| and use 
Lemma [7l D 

Remark 2. It is interesting to note that if f denotes the face with the same x- 
coordinate as f but with y-coordinate equal to n — k, the probabilities P(/ G r„) and 
P(/' G r„) being equal for symmetry reasons, we obtain that 

k n—k 

Theorem 3. As n -^ cxd, we have 

E(A„) = 4/7rlogn + o(logn) . 
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Figure 4. Zipper connection with parameter q for a fixed face /. 



Proof. Recall that ^ = ^(l+o(l)). Therefore (with a„ ~ 6„ denoting linia„/6„ = 1) 



E(A„) = J]p(/Gr„) 



/ 



n-1 



y y 



I 6 j , 6j ^ 



" X,J/=1 y^ = ly-=l 

2 n-1 J/ 2/ 



a;,y=l y,=l yj=l 
n(r2 — 1)^ f rfi n^ 



n^ 



;i+o(i)) 



2 27r 



logn + 0(ri^^ 



41ogn 



7r 



+ o(logn) 



where we have simplified (using ([6]) and some trigonometric identities) 



n-l y y 

a:^,y=i yi=i yj=i 
14 



Cj, Cj J 



n-1 y y . , , c"Kyi-\)\ ('^Kyj-\)\ 






fc=l x=l 

n-1 y ^u„, l^ J' 

X 



y=l j/i=l j/j=l 

6 "-^ ^ 2 V-^ sin2(7rA:/(2n)) + siyi^tiI / {2n)) ^8 sin 

(n- 1)2(2^-1) , ly^_ sm2(i^)/sin2(^) 



n 



2n' 



6 4 ^^^ sin2(7rA;/(2n)) + sin^ {-kI / {2n)) 

\ogn + 0(n^ 



S 2 

^ "- , ^/ 2 

loffn -I- Oir 

2 27r 



a/b , 6/a _ 1 , 1 2 



The last line can be seen as follows. Usinff that -tt. H — tt = — h r tt and 

=3 a+0 a+b a b a+b 

symmetrizing the sum in k, I we obtain that 

V sin^(i^)/sin^(|i) ^ r.-nV^ VV ' 



^^^ sin^(7rA;/(2n)) + sin^(7r//(2n)) ^ ^ ^ sin^d^) ^ ^ sin^^^) + sin 



2n' 



23 2^2 

-n logn + C(n J 

3 TT 



where the asymptotics follow from Lemmas 10 and 11 below. D 



Lemma 10. 

^^sin2(7rA;/(2n)) + sin2(7r//(2n)) " — log^ + ^(^')- 

Proof. We consider e G (0, 1). We may rewrite the sum as one fourth of the sum 
over pairs of integers {k, I) ^ (0, 0) where k and / run from —(n — 1) to n — 1. Now 
we separate the sum in two contributions. 

One for which k'^+t^ is greater than eri^. This gives a total contribution oiO{'n? /e). 
On the other part, where k"^ + /^ is smaller than en^, we may use a series expansion. 
We have 

15 



Hence, 



n— 1 n— 1 



2 



^^sm2(7rA;/n) + sm2(7r//n) vr^ ^^^^^^^ P + /2 V v ;; 

Observe the simple geometrical fact that the cardinality of the set D{r) = {{k, I) G 
Z2, A;2 + /2 < r} satisfies 7r(v^ - v^)^ < D{r) < n{^/^ + V2f. Hence /^(r + 1) - 
D{r) = TT + ©(I/a/t). Furthermore, if r — 1 < A;^ + /^ < r, we have l/(fc^ + /^) = 
l/r(l + 0(l/r)). Therefore, 



1 "^"^ 1 

E JJT? ^ E7(i + 0(-^)) 



fc2+i2<en' 



7rlog(en2) + 0(l) 
27rlogn + 0(l + loge) 



Hence, 



n—l n— 1 



^^sin2(7r/t/(2n)) + sin2(7r//(2n)) vr l^S^ + ^(^V^') ■ 

D 
Lemma 11. 

n-l ^ 2 

^^^ii^M^s^""'"^^- 

fc=l ^2?!'' 

Proof. Let p(2;) = nfc=i (-^^ ~ 4sin^(|^)). The desired quantity is —A^-Ll. Writing 

p{z) = YYk=i {^ ~ 2(1 — cos ^)) we see that p{2z + 2) has roots cos{7ik/n), that is, it 
is the Chebyshev polynomial of the second kind Ch„_i(z), which can also be defined 
by 

7 Ch„_i COS0 = — ^• 

We then have 

_,pM = g Gh;_,(-l) 
p(0) Ch„_i(-1) 

and the result follows from differentiating ([T]). D 
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4.3. Second moment of the area. Recall that 

fj' 

where the sum is over all bounded faces of Qn- 

These probabilities can be computed by Lemma |8] using the vertical zippers illus- 
trated in Figure [5j Denote by Cj the k east-oriented edges below / enumerated from 
bottom to top. Similarly, denote by e^ the k' edges below /'. When / and /' are not 
in the same column, we have 



k k' 



1=1 j=l 

When they are in the same column we rotate the square, interchanging columns and 

rows. 



f 1 

UJq' 






m. 


-^q' 


-^q 


-^q' 


-^q 


-^q' 


-^q 


-^q' 


-^q 


-^q' 


-^q 


-^q' 


-^q 


^ — 


1-^^ 



Figure 5. Zipper connection with parameters q and q' for two fixed 
faces / and /'. 



Theorem 4. As n -^ oo, 



E(A2) = Cn^ + o{n^) 



where C = 51|^ ^ .281, where 



E 



1 
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.527. 



Proof. By symmetry, the sum "^r f/^{f, f € F) over all ordered pairs of faces / 
and /' may be rewritten as the sum of three contributions: a sum over /, /' lying in 
different rows and columns; twice the sum over faces lying in same rows and different 
columns; and a sum over faces / = /'. 

Therefore, using the parametrization of faces / by the coordinates (x, y) of their 
lower left vertex, we may write 



n— 1 n—1 



x,y=l x'>x,y'>y x,y=ly=y\x'>x 



Let us compute the first sum in ([8|. Define 



, n-l n-1 y y' 

x,y=lx'>x,y'>y x,y=l x'>x,y'>y yi=l y' =1 



Recall the expression ^ for T(ej, e'j). It follows that h only depends on x, x' through 
the term 



A[n,kj:= >, sm( ) sm( ^ — 



n n 

x'>x 



by Lemma 12 So 



"-1 9,,,/ "-1 

'. = -EE^E^(".*)- 

j/=l y'>y k=l 

n-l n-l I ^■2rnk 



y,a'=i,fc,i=l \ ^2nJ \2n' yi=l yj=l 

y'>y 

18 



and using a simple trigonometric identity in the second line 



^ 12n^ -T.Mn,k) 

k=i 



sin{lTTy/n) sm{lTTy' /n 
4sin2(^) 



= ^ r^^ ^ nn-1 -2 > cot^ — 

\ k=l,kodd 

where we have used Lemma [16] in the last line. 
Now expanding the A{n,-) yields 

1 [ ^'^ sin^(f^)/sin^(|^) _ 2n V ('^^'^^^^'^^'^^0 '"^'^^^ 



A = - ^ u' y 

16^2 I ^^ 

:^ (sin2(i)/sin2(|^))cot2(i) 

n^ (sin2(|^)/sin2(^))cot2(i)cot2(|^ 
+4 > -^ 

fc,/=l,fcodd,iodd ^"M2n^^^"M2n 



^ :^ (sin^(i)/sin^(i))cot2(^) ^ -i ( 

_J_ :^ (sin^(l^)/sin^(|i))cot2(|^)cot2(^) 

4^2 Z^ siTi^r^l +sin2('^l +C(n), 

^'^^ M=i,fcodd,«odd sm ^2^j+sm ^2nJ 
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where we used that Ylki=i v"2|-llx_^^i"2c^^| = 0{n^), as shown in the proof of Theo- 



rem [3l 



sm2(f|)+sin2(|l) 



Now, using Lemmas 13, 17, and 18, we obtain that 



and hence 



^^^.n^ + Oin'), 



-n" 



a;,j/=l x'>x,y'>y 

Let us now compute the second sum of (|8|. Define 

, n— 1 71—1 y y 

^^ = ^E E n.(/./'<Er,.) = -E E EE^t^'.O 

x,y=l x'>x,y'=y x,y=l x'>x,y'=y yi=l y'.=l 

We have 

n— 1 „ 9 n— 1 , 1 I A n— 1 n— 1 n— 1 , , , 

^^ = -EEi^E-"(—)-"(—)-;^EE(E-(— )-(—)) 

x,y=lx'>x fc=l y=l k,l=l x,x'=i 

x'>x 



sin^(i 



2n{n 



,7rl 



n, —1 



1/2) , 



2cnk\ y 

J- > cos( 

^) 4sin2(|^ 



A;=l 



8 .^ sin^ri^l + sin^r^) 4n . , ^ . . sin^f?^) 






in^^lny/n) 



sin 



=in^(i)8sin2(|^) 



y^ sin^(f!)/sin^(|^) 1 y^ 

— 4^n + 0(r^^), 

-sm (2;^) 



sin2(|f) + sin2(|i) 



1 -^ (^sin^(fg)/sin^ 

2^ .. sin2(|^)+sin2(|i) 
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An ^ 

k,l=\,koAA 



where we again used that YJk,iii ''sL^{^)+sin^{^) = 0{n^). Using Lemma 
have I2 = 0{n^); multiplying by ^ we get that 

n-1 

E E P(/,rer„) = 0(n). 

x,y=l x'>x,y'=y 
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we thus 



The third term of (pi) is O(logn) by Theorem pi Therefore the only term of (pi) 
that contributes to the first-order asymptotics of the sum is the first one, and we get 
the result. D 



Lemma 12. 

E 



"-' . ifa . ifa' -n + (l-(-incot'(g: 

sm( )sm( ) = ^^ 

n n 4 

x'>x 



Proof. This is obtained through expansion in terms of exponentials. D 

Lemma 13. 

^ (sin^(i)/sin^(i)) cot^(i)cot^(i) _ 1 2^^ , 



M^Mo..o.. ^in^(i) + s-^(i) ^12 vre 



where 



k,l=l, both odd ^ 



Proof. Using that 



a^ + h'^ ^ a^ + b'^ ~ a464(a2 + 52) - ^254 + ^452 ^252(^2 + 52) 

and symmetrizing the sum in k, I, we obtain that 



^ (^sin2(ff)/sin2(^))cot2(|^)cot^ 

^ sin2(|^) + sin2(|^) 



nk. ^ cos^fl^) :^ cot2(f^)cot2(f^^ 



A:,i=l, both odd V27i^ ' V 2n 

"-1 ttA; """^ cos2(— ) 
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fc=l,feodd 2« fc=l,fcodd^^^ *^2n) fc/=l,bothodd^^^ ^2n) +^^^ (L^ 



and using Lemma 16 that 
,6 "-1 



n 
12 



E 



1 



k,l=l,hoth odd 
n-1 



■th 

1 2^/3, 



sm'(i)sin2(i)(sin2(f^) + sin2(^)) 
1 ^ 1 



E 



. 2n> 



2n' 



12 TT 



6 



sin^(i)sin^(i) ,,.,^,„,,sin^(f!) + sm^-' 



n'' + 0(n^ 



+ 0(n^ 



where we have used Lemmas 14 and 15 and the following consequence of Lemma 11 



n-l 

E 



, , , , , ,, sin2(f^)sin2(f^) 

fci=l, both odd ^2n^ "^2n' 



n-l 

E 



sin^f— ' 

fc=l,fcodd '^Sn' 



4n^ 
2" 

4^2 vr^ 



y - 



fc>l,fcodd 



n 



7r2 6 ^ 4' 
+ 0(n3). 



0(n: 



0(n, 



D 



Lemma 14. 



n-l 

E 



k,l=l, both odd 



sin^f^l + sin^^^' 



0(r2^ logn) 



>2n' 



2?!^ 



Proof. This is a consequence of Lemma 10 



Lemma 15. 

n-l 



E 



1 



6^6 



A:,J=1, 6o4/i 



^Sin2(f^)sin^(^)(sin2(if) + sin2(^)) 



2^ra' 



TT" 



E 



k,l=l, both 



,m'^{k^ + p] 



D 



+0(n^ 



Proof. The proof follows the same scheme as the beginning of the proof of Lemma [10 
Lemma 16. 



D 



n-l 

E 

k=l,k odd 



cot" 



,iTk, 



n 



2n' 2 
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+ 0{n), 



and 

Proof. The proof follows the same scheme as the beginning of the proof of Lemma [T0| 
For the first equality, we obtain that the sum is equivalent to ^C,°{2) = n^^(l — 
1/4)77^/6 = n^ /2, where C°{.s) = J2k>ikodd^~'^- -^^^ ^^^ second, we obtain that the 
sum is equivalent to ^C°(4) = n^^(l - 1/I6)7r790 = n^/6. D 

Lemma 17. 

-1 fsin2(f^)/sin2(^)')cot2(f^) 

,,4f.. sin2(f^) + sin2(fi) ^^- 

k,l=l,kodd y2n' ' ^2n' 

Proof. We use that there exists a constant c > such that ex < sin(x) < x for all 
X G [0, 7r/2] to bound the summand by (^Y i2iilj^i2\ - Since p(fe2^ p) is summable, we 
obtain a O(n^). D 

Lemma 18. 



^2n' n^ 



-i (sm^(a)/sm'(|i))cot^(£ ^„ 

Proof. Let e G (0,1) and a > 1. We decompose the sum in three contributions. The 
first one correponds to / > en. This term gives a contribution of \0(}2!kl^ sin^(^)) = 
Oin/e''). 

We now write the remaining sum (for / < en) as the sum of two contributions: 
k < al and k > al. 

This first subcontribution gives a 0{a?n'^ Ylik i>i WJW+Wii^ + ^((1 + Ci'^)^^))- 

The second gives a contribution of 

I odd, Z<£n 

k,k>al 

The first two terms are 0{n^) and the last one becomes equivalent to 
(V)'("-1) E ^ = ^(C(4)-^C(4))+0(„V) 

1=1,1 odd 

= y + 0(n^), 
since ^^l//^ = 0(l/n). D 
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5. Other domains 

We consider a planar domain (closed, simply-connected) D G C with smooth 
boundary dD. For any £ > 0, we consider a finite graph D"^ whose vertices are 
a subset oi D (1 eT?, where edges join points at distance e, and which is "simply 
connected" in the sense that the union of its closed faces is simply connected. A 
sequence (-D^) . ^ is said to approximate D if the Hausdorff distance from D^ to D is 
0{e). In this section, we give asymptotics for the two-point probabilities P(/, /' G Fg) 
and the first and second moments of the area of the uniform cycle-rooted spanning 
tree r^ on an approximating sequence (D^)£>o. The proofs are generalizations of the 
proofs of the previous section and somewhat technical so we only give sketches here. 
For complete proofs, see [5]. 

In this section, A = — J^ — J^ is the positive definite Laplacian. Recall that the 
Green's function with Neumann (resp. Dirichlet) boundary conditions is the (unique 
up to constant) smooth symmetric kernel over D solution to the PDE 

Ayu{x,y) = 6x, 

with boundary condition ^^^/V = for y G dD, where n(y) is the normal vector 
at the boundary point y (resp. u{x,y) = for y G dD). On a bounded domain 
the Neumann Green's function is not defined in general; we can only invert A^ on 
functions of mean zero. 

We denote g^j-, the Neumann Green's function and g^ the Dirichlet Green's func- 
tion. They both can be written in terms of the planar Brownian motion, respectively 
refiected or absorbed at the boundary. 

A rectilinear approximation of D is an approximating sequence (iI'^)e>o such that 
the boundaries are locally horizontal or vertical: any site in dD^ belongs to a piece 
of horizontal or vertical segment of dD"^ with length S = 6{e) such that e = o{S). We 
suppose that all our approximations are rectilinear. 

We furthermore suppose that the domain is smooth enough so that g}) is at least 
twice different iable. 

Given a point z E D and an approximation {D^)^^q, we say that a sequence of 
vertices v^ G D^ approximates z ii \z — v'^\ = 0{e). 

Let T^ be the transfer impedance on D^. 

Theorem 5. Let zi ^ Z2 be two points of D, possibly on the boundary, and v\,v\ 
approximations in D'^ of zi,Z2 respectively. Let e\,e\ be two oriented edges of D^ 
with fixed directions ei and €2, independent of e, and with starting vertices v\ and 
V2- Then 

(9) ^^«,<) =^'^^^L(^l,^2) + 0(5^), 
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where for i = 1 and 2, the symbol -^ denotes the partial directional derivative along 
direction Ci for the i-th variable. 

If we restrict to the case where Zi and Z2 are in the interior of D, then we actually 
don't need that D admit rectilinear approximations, and the proof works for any 
smooth domain. 

In order to prove Theorem [5} we actually need its analog for D = M the upper 
half-plane. We will also use this result in the proof of Theorem [7} The proof is in 
fact easier in this case, see the sketch below. 

The Green's function on C has the explicit form g^lx, y) = —^ log \x — y\- On the 
upper half-plane H, the Neumann Green's function is given by gli{x, y) = g^{x, y) + 

Now, for any infinite planar graph Q with a (possibly empty) distinguished set of 
vertices called the boundary, the Neuman Green's function is defined in terms of the 
simple random walk as follows 

oo 

Gl{x,y) = 5^(P(SRW,(n) = y) -P(SRW,(n) = x)) , 

n=0 

where SRWa;(n) is the n-th position of a simple random walk on Q, started at x and 
reflected on the boundary. 

Take Z x N to be the upper half-plane square lattice with boundary Z x {0}. Then 
we have Gz^^{x,y) = Gz'2{x,y) + Gz2{x,y). The transfer impedance is defined in 
terms of the Green's function in the same way as for finite graphs. Note that T^^ is 
the limit of the transfer impedance on a sequence of finite subgraphs exhausting Z^, 
see [1]. 

For D = C or D = M and for any e > we define D^ to be eQ where ^ = Z^ or 
^ = Z X N, respectively. 

Theorem 6. Let D = C or M.. Let zi ^ Z2 be two points of D, possibly on the 
boundary, and fi,ff approximations in D^ of zi,Z2 respectively. Let el, 62 be two 
oriented edges of D^ with fixed directions ei and 62, independent of e, and with 
starting vertices vf and ff . Then 

(10) T%el^,elJ=e'^^^^g^^{z,,Z2) + o{e'), 

where for i = 1 and 2, the symbol -^ denotes the partial directional derivative along 
direction Cj for the i-th variable. 

25 



Sketch of proof of Theorem \^ We start by showing ( 10 ) in the case where D is the 
whole plane using the asymptotics 

Gz2(0, z) = -— log \z\+c + 0{-^) , 
2n \z'^\ 

for the Green's function on Z^, see [TT] , knowledge of the error term (which extends to 
a continuous harmonic function), and the mean value property estimates (Proposition 
A. 2 of ^) to control the remainder terms. 

We then deduce the formula for the half-plane by using a reflection principle, that 
is, writing that 

G'zxn(«,^) =Gj2{u,v) + Gz2iu,v) , 
where G^xn is the Green's function with Neumann boundary condition on the half- 
plane square lattice. D 

Sketch of proof of Theorem^ Recall the definition T{ab,uv) = G{a,u) — G{b,u) — 
G{a, v)+G{b, v) of the transfer impedance in terms of the Neumann Green's function. 
The proof follows by a careful study of the error term when replacing the discrete 
Green's function by (?£,. In the case where both zi and Z2 are interior points, we 
may use the convergence results of discrete Green's functions to their continuous 
counterparts and the mean value property estimates ^ as follows: 

Let D' be a compact simply connected subset of the interior of D containing both 
Zi and Z2. For any Wi in D', define 

where uf is the oriented edge with direction ei whose starting vertex is wi. 

We now use the mean value property estimates for the variable wi to replace each 
term of H"^ by its average over a discrete ball of small radius r > 0. For fixed wi and 
r, the limit lime_>o H^{wi) exists and is equal to 



— 2(1 ■^—hD{wi,Z2)ds~ — - 



-hDiwi, Z2)ds ~ j—hoiwi, Z2)ds] , 

ldB{wi,r)+ *-'^2 JdB{wi,r)~ ^^^2 ' 

when w\ G dD' , where hY){w\,W2) = g''j^{wi,W2) — gc{'Wi,W2) and dB{wi,r)~^ and 
dB{wi, r)~ are the right and left half-circles the union of which is the circle of radius 
r around Wi. 

By further taking r — )■ 0, we obtain that lim£_^o -^^^(""^1) = gfge ^£'('"^1; -^2) when 
wi e dD'. 

We now let g be the harmonic extension in D' of the function gg ge ^d{wi^ Z2) on 
dD' . It follows from classical results that for any wi G D' we have H^{wi,Z2) — 

6, we 



g{wi,Z2) = 0(1). Since ^T£Z2(M^,e|) converges to g^^gc{wi, Z2) by Theorem 



then deduce that ^T^(Mf , e|) converges, as a function of wi, to a function which takes 
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boundary values g^ 8e dhi'^i^ ^2) and has a single second order pole with coefficient 
— ^. This uniquely determines its limit to be gfg^ g}){wi, Z2) on the whole domain 
D' . In particular, we have 

1 92 

In the case where at least one of the points Zi or Z2 is on the boundary, we use a 
reflection argument (making full use of the rectilinearity of the domain), and use the 
result for interior points. D 

In a way, the previous theorem is a refinement of a result of [U] (Lemma 17), where 
the quantity -{G{vl + £,^2) — Giyl.vD) was shown to converge to the ffist partial 
derivative of (7^. 

We now suppose that D is bounded. Let n^ be the number of spanning trees of 
D^ , and Ag its number of cycle-rooted spanning trees. 



Theorem 7. Let zi and Z2 be two distinct points in D, possibly on the boundary, 
and vl.vl G -D^ sequences of vertices approximating zi, Z2 respectively. Let f^e and 
/^e be the faces whose lower left vertices are ff and ff . Then 

F(/.j,/.|Gr,) = ^(^°,(^i,^2) + o(l)) 

Sketch of proof. We do the proof in two steps. 

First, we use Lemma |8] to write this probability as a sum over a product of two 
paths of edges. We let 7^ and 72 be two rectilinear, non-crossing dual paths in 
D'^ from the unbounded face to /j,e and /^e, chosen in such a way that these paths 
approximate (with 0{e) error) two non-crossing continuous paths 71 and 72 from dD 
to Zi and Z2- Lemma |8] states that 



K 

eie£(7f),e26£(7|) 



iei,e2j 



where -E(7i) and £^(7!) are the set of coherently oriented edges (let us say, counter- 
clockwise around /„e and f^^^) on the paths 7f and 7|. Using Theorem [s] we thus 
obtain, by a convergence argument, that 

(11) -P(/.^/.i) = - / / ^^gUwi,W2)\dwr\ \dW2\+0{l). 

where for i = 1,2, notation e, denotes the normal direction to the right of path 7^ 
when followed from dD to Zi. 

The second step is to evaluate this integral. We observe that this integral trans- 
forms covariantly under conformal transformation: if : D — )■ EI is a Riemann map 
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to the upper half-plane EI then 

1 



gniwi,W2) 



27r 



log 



(Wi) -(J){W2)\ - ^ log I 



(Wi) -(I){W2)\. 



This is the real part of the function 

goiwi, 102) = -^;^ log(0(wi) - (t){w2)) 



1 



27r °'"' ' "' '' 27r 
Moving the integration to the upper half-plane with z. 



log(0(wi) -0(^2))- 



[Wi 



Re 

27r 



27r 



log 



H11) -'Hii) 



Zi - Z2 



_d__d_ 

dzi dz2 



\og{zi - Z2)dZidZ2 - 



(11) becomes 
d d 



</)(7l) J(j}{j2) 



dzi dz2 



Z\ - Z2 



where we used the conformal invariance of g\) in the last equality. 
Let Lg and A^ be the combinatorial length and area of V^. 



log(zi — Z2)dz\dz', 



D 



Theorem 8. We have 



and 



EfL, 



EfA, 



= 8 + 0(1), 
loge + o (loge) 



^{Ki)=e-'C{D)\D\\o{e-') 



where \D\ is the Lebesgue measure of D and C{D) = ^§12 fjj2 gDi^^'^)\d'Z\'^\dw\'^ is 
the mean normalized exit time referred to in the Introduction. 

Sketch of proof. The first equality is an easy consequence of a result of Levine and 
Peres [I2]. The two others follow from the fact that moments are obtained as sums 
of probabilities as given in Lemma [3j 

In order to evaluate the sums, we use that 



T'ie'eV 



-15 ^2) 



^eZxeN(ei,e2; 



< Cos'' , 



for some Co > 0, which we show by the mean value property estimates and a reflection 
argument near the boundary. This allows us to replace T^ by T^zxen with an error 
term of O(e^). 

The second equality follows by using this approximation and an explicit compu- 
tation for T^zxsN- 

For the third equality, we show convergence of the sum "^f f> IP(/, /' G T^) to an 
integral. This works without trouble for the contribution of pairs of faces that are 
macroscopically far apart and follows from Theorem [7j For the pairs of faces that 
are within distance 0{e) of one another, we replace T^ by T^zxeN if they are in the 
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bulk, and if they lie on the edge, we use the electric flow interpretation to show that 
their contribution is negligible to the first order. D 

The expected length and area do not depend on the shape of the domain. However, 
the mean normalized exit time (and thus the second moment of the area) does depend 
on this shape. Here are two examples. 

(1) Let D be a 1 X r rectangle. We have C{D) = 5i?^^ where 



/3(-) = E 



k,l=l, both odd ^ ^ 

This follows from Theorem [8] using the explicit form of the Dirichlet Green's 
function on Qn,nT, the n x nr rectangle in the square lattice (computation 
omitted). 
(2) For the unit disk C{D) = l/vr: the Dirichlet Green's function is —^ log | fE^ \ ■ 
Substituting z = rie*^^ and w = r2e*^^ we get (using symmetry we assume 
^2 = and multiply the result by 27r) 



C = -TT^ log 



\D\ ./o Jo Jo 



2tt 



Tie ^ — r2 



rir2 dOi dri dr2 



- rir2e*^i 
and writing the log as the log of the numerator minus the log of the denomi- 



nator gives 

16 

TT 



I I (logmaxjri, r2} — log l)rir2 dri dr2 = —■ 

Jo Jo TT 



Polya proved that C{D) (in a slightly different form) is maximized for the disk 
(and only the disk) and the argument goes as follows: 
Let -pTjy: = inf J w{f) be the infimum, over all smooth functions / over D vanishing 

on the boundary, of w{f) := -^ — -2- This defines a quantity P{D), which Polya 

4(/d /) 

called the torsional rigidity of the cross-section D, and which is, in mechanical terms, 
a measure of the resistance to torsion of a cylindrical beam with cross-section D. In 
fact the infimum oiw{f) is realized and we have -pky- = w{fo), where /o is the solution 

of A/ = 2 with Dirichlet boundary conditions. Since fo{z) = 2 f^g'l){z,w)\dw\'^, we 
have 

\^\ Jd^ \^\ Jd 

4 {Jr>fo{w)\dw\'f _ 4 [jj,fo{w)\dw\^f _ 2 p^^^ 
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In [H], Polya proved an old conjecture of Saint- Venant (1856) that for fixed area, 
the disk is the domain that maximises P{D). Therefore, since C{D) is easily shown 
to be invariant under dilation, we have proved that C{D) is maximized for the disk. 

For the benefit of the reader here is a sketch of Polya's proof which goes in two 
steps: He first shows that C{D) < C{D') where D' is any Steiner symmetrization of 
D. (A Steiner symmetrization of D along line L is a set D' having L as symmetry 
axis and such that for any line V perpendicular to L the set Lr\D' is an interval of 
length equal to the length of the set L' fl D.) He then uses a result of Gross (1917) 
(refined by Lyusternik) which states that for any domain D there is a countable 
sequence of Steiner symmetrizations whose limit is the disk with the same area as 
D. 

6. Higher moments of the area 

For any sequence {Qn) such that -Qn approximates a domain D, we may write 
higher moments of the area of Tg^ as 

E(A^j= Y. P(/i,---,/.Grgj, 

/ir-- Jk 

where the sum is over all ordered fc-tuples of bounded faces of Qn- 

For any k fixed distinct faces /i, ■ ■ ■ , fk, "we may try to compute the probabilities 

P(/i!'"" ) /fc £ ^Gn) by setting zippers from each of the fi to the boundary (or 

equivalently, branch-cuts for the Green's function) with parallel transport g^ (or 

equivalent ly, monodromy qi around faces fi). 
We conjecture the following. 

Conjecture 1. For any integer k > 2, there exists ak = ak{D) > such that, when 
n — )■ oo, 

KiAl) = akn''-\l + 0(1)) . 

7. Application: size of fjords in the uniform spanning tree 

The union of a spanning tree of a graph with any edge in its complement defines 
a unique cycle, which we call a fjord of that tree. We use the previous techniques to 
give a lower bound on the first and third moments of the length of a random fjord 
in a uniform spanning tree. 

Let Qn C Z^ be a sequence of graphs such that -Qn rectilinearly approximates a 
bounded domain D G C Consider the uniform spanning tree (UST) on Qn- Add 
an independent random edge to the tree. Then one obtains a CRST with weight 
proportional to the length of its cycle. Using the previous results, we can obtain 
information on the size of this cycle. We denote by T„ the random CRST thus 
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obtained, by ]Eust+ the corresponding expectation, and by L-f„ and Ax„ the length 
and area of its cycle. 

Lemma 19. For any measurable function f , we have 

E(/(r„)L„) 



EusT+(/(T„)) 



EfL, 



-'nj 

Therefore, the results on the moments of the area of the uniform CRST proved 
earlier show that 



Theorem 9. 



Iqct ■ri 

Eust+(At„L^iJ = ^ + o(logn) 



C{D)\D\^^, , ^^,^ 
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and 

Eust+(A2^„L^;^) = ^^^g^'^' n^ + o{n\ , 

where C{D) is the mean normalized exit time. 

The isoperimetric inequality in the square grid reads L^ > 16A. As a corollary, 
we thus obtain bounds on the first and third moments of L-f . 

Corollary 1. 

o 

Eust+(Lt„) > -logn + o(logn) 

n 

and 

Eust+(LtJ > 32C(D)|D|n2 + o{n^). 

8. Questions 

1. Can one do the previous computations on other lattices? 

2. Here is a way to sample a CRST: suppose the edges of the graph are or- 
dered; sample a uniform spanning tree (using David Wilson's very efficient 
algorithm) and add a random edge from the complement. This creates a 
cycle-rooted tree. If the edge added has minimal index in the set of edges 
belonging to the cycle thus closed, then output this cycle-rooted tree. Oth- 
erwise, throw this configuration away and sample again. This method works 
but is not very efficient. Is there an efficient way to sample uniform cycle- 
rooted trees? 

3. What can be said about the distribution of the cycle of the uniform CRST 
on Z^ conditioned on surrounding the face at the origin? 
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